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bJQ. Abstract 

We argue that the conjectured dark mater in the Universe may 
^ ' be endowed with a new kind of gravitational charge that couples to 

a short range gravitational interaction mediated by a massive vector 
field. A model is constructed that assimilates this concept into ideas of 
current inflationary cosmology. The model is also consistent with the 
observed behaviour of galactic rotation curves according to Newtonian 
dynamics. The essential idea is that stars composed of ordinary (as 
opposed to dark matter) experience Newtonian forces due to the pres- 
ence of an all pervading background of massive gravitationally charged 
cold dark matter. The novel gravitational interactions are predicted to 
have a significant influence on pre-inflationary cosmology. The precise 
details depend on the nature of a gravitational Proca interaction and 
the description of matter. A gravitational Proca field configuration 
that gives rise to attractive forces between dark matter charges of like 
polarity exhibits homogeneous isotropic eternal cosmologies that are 
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free of cosmological curvature singularities thus eliminating the hori- 
zon problem associated with the standard big-bang scenario. Such 
solutions do however admit dense hot pre-inflationary epochs each 
with a characteristic scale factor that may be correlated with the dark 
matter density in the current era of expansion. The model is based on 
a theory in which a modification of Einsteinian gravity at very short 
distances can be expressed in terms of the gradient of the Einstein 
metric and the torsion of a non-Riemannian connection on the bun- 
dle of linear frames over spacetime. Indeed we demonstrate that the 
genesis of the model resides in a remarkable simplification that occurs 
when one analyses the variational equations associated with a broad 
class of non-Riemannian actions. 



1. Introduction 

The standard cosmological paradigm offers a fertile domain for the fusion of 
ideas from astrophysics and particle physics. Certain prejudices inherent in 
the classical Priedmann cosmologies are thought to find a more natural res- 
olution in terms of dynamical consequences of phase transitions in the early 
Universe leading to a variety of inflationary scenarios. Among the claims of 
success are the resolution of ad hoc fine-tuning conditions and the puzzle of 
the observed homogeneity and isotropy of matter. However a further prob- 
lem is the strong indication that the observed matter may be but a small 
fraction of the gravitating matter in the Universe. In particular without a sig- 
nificant contribution from directly unobserved matter it is hard to reconcile 
Einstein's theory of gravitation (and its Newtonian limit) with the dynam- 
ics of the Universe. Although knowledge of conditions in the early Universe 
by their very nature require huge extrapolations from the observed data, the 
conjecture that hidden matter also exists seems necessary to explain the New- 
tonian dynamics of both individual galaxies and galactic clusters. For a star 
in a circular orbit of radius r with tangential speed v within a spherical dis- 
tribution of galactic matter with mass M(r), Newtonian dynamics predicts 
terms of the Newtonian gravitational coupling G. Suppose 
Mq = M(ro) is the constant mass of the visible galaxy where Tq denotes the 
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radius of visible matter. If r > ro then in the absence of gravitating matter 
outside the visible galaxy v^{r) = i.e. v{r) ~ This however is 
not observed in general. For r > ro, indications are (e.g from Ha emission) 
that f ^(r) = k for some constant k, or M(r) ~ r requiring a matter density 
distribution p(r) of the form p(r) ~ The material in M(r) for r > Tq has 
been called dark matter. 

There have been many suggestions for the origin of these observations 
including a number that contemplate modifications of Newton's law of grav- 
ity and the effects of new interactions ||, 0, ||, ^ ^, §]. This paper 
explores the idea that they result from the existence of dark matter that is 
characterised by a novel gravitational interaction and that the carrier of this 
component of gravitation is a vector field that influences the pre-infiationary 
phase of the Universe (c.f. As such it may modify the inflationary period 
needed to address the so called "flatness problem" and may put the so called 
"horizon problem" into a new perspective. Depending on the sign adopted 
for the coupling of the Proca field in the fundamental action defing the theory 
this interaction may also be responsible for removing the classical curvature 



singularity associated with the traditional Friedmann cosmologies [|T0|- The 
coupling of the dark matter to the Proca component of gravitation is defined 
by a new kind of conserved charge, the value of which can be correlated with 
the "mass" of the Proca quantum and the behaviour of the scale factor of 
the Universe. 

An important feature of this description is that it is rooted in recent de- 
velopments in non-Riemannian descriptions of gravitation. Non-Riemannian 
geometries feature in a number of theoretical descriptions of the interac- 
tions between fields and gravitation. Since the early pioneering work by 
Weyl, Cartan, Schroedinger, Trautman and others such geometries have of- 
ten provided a succinct and elegant guide towards the search for unification 



of the forces of nature |]rT| . In recent times interactions with supergravity 
have been encoded into torsion fields induced by spinors and dilatonic in- 
teractions from low energy effective string theories have been encoded into 
connections that are not metric-compatible |jl2|, |T3|, 0. However theories 
in which the non-Riemannian geometrical fields are dynamical in the absence 



of matter are more elusive to interpret. We have suggested elsewhere [|16 
that they may play an important role in certain astrophysical contexts. Part 
of the difficulty in interpreting such fields is that there is little experimental 
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guidance available for the construction of a viable theory that can compete 
effectively with general relativity in domains that are currently accessible to 
observation. In such circumstances one must be guided by the classical solu- 
tions admitted by theoretical models that admit dynamical non-Riemannian 

This approach is being currently 



structures O, K m 20, 21, 22, 23, 24 



pursued by a number of groups |jT6|, ^ 



The 



work below continues in this vein by recognising that the essence of the non- 
Riemannian approach is the existence of new gravitational interactions. The 
nature of a particular interaction is explored in the context of modern infla- 
tionary cosmologies and the dark matter problem. Although the dark matter 
model below is formulated in terms of the standard Levi-Civita Riemannian 
connection the introduction of the Proca component of gravitation is best 
appreciated in terms of its genesis within non-Riemannian geometry. 



2. Non-Riemannian Geometry 

Einstein's theory of gravity has an elegant formulation in terms of 
(pseudo-) Riemannian geometry. The field equations follow as a local ex- 
tremum of an action integral under metric variations. In the absence of 
matter the integrand of this action is simply the curvature scalar associated 
with the curvature of the Levi-Civita connection times the (pseudo-) Rie- 
mannian volume form of spacetime. Such a connection V is torsion-free and 
metric compatible. Thus for all vector fields X, Y on the spacetime manifold, 
the tensors given by: 

T(X, Y) = V,Y - V,X - [X, Y] (2.1) 

S = Vg (2.2) 

are zero, where g denotes the metric tensor, T the 2-1 torsion tensor and S 
the gradient tensor of g with respect to V. Such a Levi-Civita connection 
provides a useful reference connection since it depends entirely on the metric 

o 

structure of the manifold. This special connection will be denoted V in the 
following. 

A general linear connection V on a manifold provides a covariant way 
to differentiate tensor fields. It provides a type preserving derivation on the 
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algebra of tensor fields that commutes with contractions. Given an arbitrary 
local basis of vector fields {Xa} the most general linear connection is specified 
locally by a set of 1-forms A"^ where n is the dimension of the manifold: 

V,^Xt, = AMXa)X,. (2.3) 

Such a connection can be fixed by specifying a (2, 0) symmetric metric tensor 
g, a (2-antisymmetric, 1) tensor T and a (3, 0) tensor S, symmetric in its 
last two arguments. If we require that T be the torsion of V and S be the 
gradient of g then it is straightforward to determine the connection in terms 
of these tensors. Indeed since V is defined to commute with contractions and 
reduce to differentiation on scalars it follows from the relation 

X{g{Y, Z)) = S{X, Y, Z) + g(V,y, Z) + g{Y, V,Z) (2.4) 

that 

2g(Z,V,r) = X{giY,Z)) + Y{g{Z,X))-Z{g{X,Y)) 

- g(X, [Y, Z]) - g{Y, [X, Z]) - g{Z, [Y, X]) 

- g{X, T{Y, Z)) - g{Y, T(X, Z)) - g(Z, T(F, X)) 

- S(X, F, Z) - S(F, Z, X) + S(Z, X, Y) (2.5) 

where X, F, Z are any vector fields. The general curvature operator Rx,y 
defined in terms of V by 

^x,yZ = V,V,Z - V,V,Z - Vj,,,,Z (2.6) 

is also a type- preserving tensor derivation on the algebra of tensor fields. The 
general (3, 1) curvature tensor R of V is defined by 

R(X,y,Z,/3)=/3(i?,,^Z) (2.7) 

where j3 is an arbitrary 1-form. This tensor gives rise to a set of local curva- 
ture 2-forms R""},: 

i?%(X,F) = ^R(X,y,Xfe,e'^) (2.8) 

where {e*^} is any local basis of 1-forms dual to {X^}. In terms of the con- 
traction operator ix with respect to X one has ix^ = 1^6" = e°'{Xi,) = 5";,. 
In terms of the connection forms 

R\ = dA\ + A\ A A\ (2.9) 



5 



In a similar manner the torsion tensor gives rise to a set of local torsion 
2- forms T": 

r"(X, y) = ^ e''(T(X, Y)) (2.10) 

which can be expressed in terms of the connection forms as 

^« = de" + A^Ae^ (2.11) 

Since the metric is symmetric the tensor S can be used to define a set of local 
non-metricity 1-forms Qab symmetric in their indices: 

Qab{Z) = S{Z,Xa,Xb). (2.12) 



3. Non-Riemannian Actions 

To exploit the geometrical notions above in a gravitational context one 
must construct field equations that determine the torsion and metric gradi- 
ents. These arc most naturally derived from an action principle in which the 
metric, components of the connection and matter fields are the configuration 
variables. One requires that an action be stationary with respect to suitable 
variations of such variables. If the action 4-form yl(g, V, ■ • •) in spacetime 
contains the Einstein-Hilbert form 

AEHig,v) = n^i (3.1) 

where * denotes the Hodge map on forms and Aeh denotes the variational 

g 

derivative of Aeh{s, V) with respect to g, then 

Aeh = -h""' Ein(X„, X,) * 1 (3.2) 
g 

where hab = g{Xa,Xb) and 

Ein = Ric-^g7^ (3.3) 
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is given in terms of the Ricci tensor Ric(X, Y) = R(Xa, X, Y, e") by 

mc{X, y) = ^ (Ric(X, Y) + Ric(F, X)) (3.4) 

and TZ = Ric(Xa,X"). Further details of these variational calculations may 
be found in |T^. Unlike the Einstein tensor associated with the Levi-Civita 
connection, the Einstein tensor Ein defined in ( p.3|) ) is associated with the 
non-Riemannian connection V and is not in general divergenceless: 

(V ■ Ein)(Xfe) = (Vx.Ein)(X",X6) ^ 0. (3.5) 
It is therefore instructive to decompose the non-Riemannian Einstein tensor 

o 

Ein into parts that depend on the Levi-Civita connection V. For this purpose 
introduce the tensor A by 

A(X, Y, = (3{VxY) - (3{VxY) (3.6) 

for arbitrary vector fields X, Y and 1-form (]. In terms of the exterior covari- 
ant derivative |]31| Dand the Ricci tensor Ric associated with the Levi-Civita 



connection, one may write: 

Ric(X„ Xfe) = R°ic(X„ X,) + i„ i, (DA^b + A^^ a (3.7) 

where A^^ = A(— , X^, e") is a set of local 1-forms. In terms of these forms 

T'' = X\r.e' (3.8) 

Qab = — Aafe — Afea- (3.9) 

It follows from ( |3.6| ) that Ein differs from the Levi-Civita Einstein tensor 

o 

Ein by terms involving the tensor A and its derivatives. However for a large 
class of actions containing the torsion and metric gradient fields one finds 
that these terms can be dramatically simplified. 

To see this simplification most easily it is preferably to change vari- 
ables from g, V to g, A in the total action 4-form and write A{g, V, ■ ■ ■) = 
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£(g, A, ■ ■ ■). Since V depends only on the metric it follows from (|3.6| ) that 

o 

^_A_^ = V and ^A^ = — ^_V^. Hence the variational field equations are 

^ ^ ^ ^ = ^ = (3-10) 

V A 

^ = ^ + ^^ = 0. (3.11) 
g s s X 



From ( p.lOp one sees that the term C in (|3.11|) does not contribute. To 

A 

evaluate these variations one first expresses the action in terms of g and 
the components of A and its derivatives. To this end it is convenient to 
introduce the (traceless) 1-forms A'^f, = A'';, — ^ 5°'b A'^d, and the (traceless) 

0-forms A'^bc = ^^hc — \ S°'b ^'^dc where A^^c = ic A'^;,. If we express the total 
action 4-form C{g, A, ■ ■ ■) as 

C{g, A, ■ ■ ■) = CEHig, A) + C,{g, A, • ■ ■) (3.12) 
for some form Ci{g, A, ■ ■ ■) where 

Ceh = n- X\ A A^b A ^ (e' A ea) - d(\\ ^ (e\ e J) (3.13) 

o 

in terms of the Levi-Civita scalar curvature TZ then 

Ceh = -h""' Ein(X„ X,) ^ 1 - h^' Sab (mod d) (3.14) 



where 

Sab = ^ A^d A A'^p A {5'afe * (c^ A Cg) - 6^ a * (Cfe A Cg) - ^ (Ca A Cg) } . (3.15) 

Thus the Einstein field equation is 

Ein(X„, Xfe) * 1 + Sab + Tab = 
where Ci = —h"-^ Tab- Next the variations with respect to A yield 
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Ceh^ = A {A'^a A * (e\ -X^ca * (e'' a Ca)} (mod d) 
A 

\ a 

''-1 — ^ fe A a- 

A 

By splitting off the trace part of the resulting field equation one may write: 

= (3.16) 

*(e\e,)-?eA *(e'=Aea)+^''a = (3.17) 

where J-'°'b = ^""b — \S"'b^'^c- To illustrate how the field equations (|3.16|) and 
( |3.17| ) can greatly simplify the terms in ( p.l5| ) consider the eight parameter 



class of models in which the torsion and metric gradient fields enter the action 
according to: 

Ci=4KR\^^R\-2i^l + aiQ;,^Q 

where k, ak are arbitrary coupling constants and ^ is a cosmological constant. 
The 1-forms u and v may be expressed in terms of the torsion forms T*^ and 
non-metricity forms Qab as follows: 

v = \a\e'' = \Q-\Q-T 

where T = i^T", Q = = -2 A^^ and Q = e'^i^'Qab- Furthermore = 
— |d(5 is proportional to the Weyl field 2-form dQ. Computing the variational 
derivatives above one finds that ( |3.16| ) yields: 

d^dQH — -kQ-l — -ku-\ — -kv = 

lOK lOK lOK 

(3.19) 

while ( p.lTp implies: 

u = (3iQ (3.20) 
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and 



v = P2Q (3.21) 

A be — ^^j C 6 Ic V + 

12 



24 



(3.22) 



^ 6 as + 3 0^4 — 21 tts — 3 ae + 54 0304 — 27 as^e — 2 
16 — 6 q;2 — 6 03 + 42 ag + 27 — 108 0203 

^ 6 02 — 21 ^4 + 3 cts — 3 tte + 54 a;2a;5 — 27 a^a^ — 2 

16 — 6 ^2 — 6 03 + 42 ag + 27 a^^ — 108 a;2a3 

Substituting (|3.20|) and ( p.21| ) into (|3.19|) one now obtains a Proca-type equa- 
tion in the form 

d^dQ + /33*(5 = (3.23) 

where 

16q!i — q;4 — as 16a4 — 2a2 — IQa^ — 2a3 — oq 

^3 = 77 ^ 77 Pi 77 P2. 

16/? 16k 16fi: 

From the metric variation of Ci one finds 

Tab = igab-^'^ + T{Xa, X^) -k I + fab 

where 

1 1 
T = -k-^iPs {iaQ A -^hQ-'Qab Q A *Q)+iad(5 a *ifedQ-- 5(„6 dQ A ^dQ} e"(g)e* 



and 



^fc — Cl 9abQ A * Q + C2 iaQ A * IfoQ 

in terms of 



^ 13A , , 29A/^2 23A^ _ , ^ /o 

324 324 162 324 324 1296 ^ ' ^ 
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and 



2/^2 37/31^ 11 /^s' 1 

C2 = H • (3.25) 

^ 81 81 162 81 162 648 ^ ' 

Remarkably Tab exactly cancels £ab leaving the Einstein equation in the form 

Ein + £g + KT = 0. (3.26) 

Thus the solutions to the field equations derived from the action 
-^-E/flg, A) + JF where JF is given by (|3.18| ) may be generated from solu- 
tions to the Levi-Civita Einstein-Proca system (|3.23| ), ( p.26| ) with arbitrary 
cosmological constant ^. Once Q is determined the full connection follows 
from 

\\ = \\-l5tQ (3.27) 



where X^b is given by ( p.22 ). The torsion and metric gradient then follow 
from (B.8|) and 



The above discussion has concentrated on the gravitational sector of the 
theory. The inclusion of matter into JF is straightforward in principle. For 
example if the matter is minimally coupled to the Weyl form Q, the action 
T will contain a term Q a j where the 3-form j is the Proca charged matter 
current. Then ( |3.23| ) will contain an extra source term proportional to j while 
( |3.26|) will contain an addition stress tensor from the variation of the matter 



action with respect to the metric tensor. 



4. Divergence Conditions 

The previous section dealt with a general framework for deriving the vari- 
ational equations for matter interacting with non-Riemannian gravitation. In 
the following we construct models in which the new component of gravitation 
is described in terms of a Proca vector field with the matter fields considered 
as thermodynamic fiuids. Action principles for relativistic fiuid models in- 
teracting with gravitation require the imposition of sufficient constraints to 
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maintain consistency. Fortunately it is possible to write down the coupled 
field equations for an Einstein-Proca-fluid system without much effort by ex- 
ploiting the Bianchi identities to ensure that the fluids satisfy the required 
relativistic continuum field equations. This procedure relies on recognising 
how the metric variations of the components of the action are correlated with 
the matter variations. 

Suppose the underlying theory of gravitational and matter fields (with 
compact support) is defined by an action 4-form A on spacetime M as a sum 
of 4-forms J2r-^r- Then in general there exists a set of symmetric second 
degree tensors {%}, each member of which is constructed from an element 

o 

As belonging to a subset {A^} contained in {A^}, satisfying V- 7^ = in 

o 

terms of the Levi-Civita connection V. The subset {A^} is defined so that 
each member A^ depends on a set of tensor fields {^^k} (and possibly their 
derivatives) where the intersection of the sets contains only the metric 

tensor on M. 

The essence of this result is well known (although perhaps not widely 
appreciated) and relies on the behaviour of the action functional S in response 
to variations of fields induced by local diffeomorphisms of the manifold M. 
For any set of 4-forms {A^} the action 

5= / A = ^/ A, 

where A is constructed from the metric tensor g and tensor fields of arbitrary 
type {$1, $2, ■ ■ •} on M. For any local diffeomorphism on M generated by a 
vector field V 

I £vAr= [ d(ivK)^ [ iyA^-0 (4.1) 

JM JM JdM 

since the Lie derivative £v = ^\/d + (\.iv in terms of the contraction operator 
iy, dAr = and all fields in A^ are assumed to have compact support. To 
explicitly construct the set {7^} write the variations of A^ with respect to 
e {^'k} as 

Ar — ^k ■ J'k (mod d) 

for some tensor Tk contracted with ^k- In an arbitrary local frame {Xq} 
with dual coframe {e^} such that e^(Xc) = 5^c, the symmetric "stress" tensor, 
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% = %ab^"' ® e*, associated with is obtained from the metric variation: 

= g{Xa, Xb) * 1 (mod d) . 

g 

Since the Lie derivative is a derivation on tensors, £v^r may be expressed 
in terms of the field and metric variations: 

+ {£v^i) ■J^i + {£v^2) ■J^2 + --- (modd) 

If the set {$1, $2i • • •} is unique to the 4- form then the equations {J-^ — 0} 
constitute a set of variational field equations. If these are satisfied then 

£vAr = (£vg)(X,,X,)%''' * 1 (mod d). 

Since 

i£vg)iXa,Xh)=iaVx,V + hVx^V 

where V = g{V, —) and g is symmetric one may write 



1 f £vAr^ [ r/^i„Vx,1^*l. 

2 Jm Jm 



For any % it is convenient to introduce the associated 3-form J^y by 



Now 



d{Jrv) ^ ^{TrabV^ * 6^) 

= ( Vx, %) {Xa, X„) y« g'^ ^1 + T/^ i, Vx. V <5^„ ^ 1 

^{v-%){v) ★i + r/^iA^i^*! 
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where D above denotes the exterior covariant derivative associated with V. 
Thus for fields with compact support 

I I £yA= f r/sA,t/^i 

= f dJv- f (V-T)(\/) *1 

Jm Jm 

= I JV- I {^-TrW) 

JdM Jm 

Since Jy = on dM and V is arbitrary it follows that 

V- Tr = 0. 

This result follows from the definition of % and the imposition of the field 
equations and should not be confused with a "conservation law" in an ar- 
bitrary spacetime. In general, "conservation laws" owe their existence to 
further conditions. For example if there exists a vector field V = K where K 
satisfies the Killing condition {CkS = 0) then both terms on the right hand 
side of ( [4 .21) are zero and JrK defines a genuine conserved current: 

dJrK = 0. 

Thus we have shown that every component 4-form A,, that gives rise, by 
variation, to a set of field equations {J-'k = 0} for any of the dynamical tensor 
fields in the action, excluding the metric g, also gives rise to a divergenceless 
second rank tensor %. 



5. The Einstein-Proca-Matter System 

Based on the reduction of the non-Riemannian action to a theory of gravity 
in terms of the standard Levi-Civita torsion free, metric compatible connec- 

o 

tion V, we now construct a model of gravity and matter that includes the 
Proca field in the gravitational sector. As befits its origin in terms of purely 
geometrical concepts the Proca field is regarded as a gravitational vector field 
that is expected to modify the gravitational effects produced by the tensor 
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nature of Einsteinian gravity. The model is based on ( p.26| ) in section 3 with 
i = and the right hand side replaced by the stress tensors for fluid matter. 
We may normalise the Proca field Q so that the term k t becomes aS where 
E is the Proca stress tensor, a = ±1 and the constant (3^ is replaced by the 
square of the (real) Proca field mass rria. In the following we explore the 
consequences of both signs for a. With a = —1 and = the model is 
analogous to the classical Einstein-Maxwell- fluid system. With a = -|-1 an 
interpretation of a S as a stress tensor for matter would be "unphysical" . In 
the absence of other matter and in a weak field limit it might be expected to 
lead to stability problems analogous to fermionic stresses in the absence of 
quantisation. However whatever criteria are needed to maintain the stability 
of the Minkowski background in the presence of gravitational (metric) waves 
might be similarly invoked in the presence of gravitational Proca fields. We 
shall not immediately discard the a = +1 possibility since it has interesting 
consequences for cosmology. Henceforth we restore the physical constants 
G, c, h and adopt MKS units when comparing with observation. 

In general we consider matter to be composed of ordinary matter defined 
to have zero coupling to the Proca field and "dark matter" defined to have a 
non-zero coupling. We call this coupling Proca charge and denote the basic 
unit of Proca charge by q. For our discussion of cosmology we model both 
types of matter by fluids with standard stress tensors % and Tq respectively. 

o 

Denoting the standard Levi-Civita Einstein tensor by Ein and the contri- 
bution of the Proca field a to the Einstein equation by cxS, cr = ±1 we 
have 

Ein + aS = ^(ro + r,) (5.1) 

where 

% = (cVo + Po) Vo0Vo + Pog 

Tg={C^Pg + Pg)Vg®Vg + Pgg 

g(K),K)) = -l 
g(K,K) = -l. 
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The mass densities pf, / = 0, g of ordinary and dark matter respectively are 
functions of the particle densities n/ and the entropies s/ per particle: 

Po = Po{no, So) 

Pq = Pqinq,Sq). 

The pressure Pf and temperature Tf may be derived from Gibb's relation 

c'^dpf = c^ Hfdnf + UfTfdsf (5.2) 

where Pj = ^^^^ is the associated chemical potential. In terms of the 
notation introduced in Section 3 we have rescaled the Weyl 1-form Q to the 
1-form a and r to 

^ = (a ® « - ^ a{a) ^ + (\^F ® F {F ^ ^ F) g) 

where F = da is the Proca field strength. Equation ( |5.1D for the metric 
must be supplemented by field equations for the Proca field a and the fluid 
variables together with their equations of state. In view of the comments in 
Section 4 we adopt as matter field equations 

V- To = (5.3) 



V.(^r,-.E)=0 (5.4) 

o 

which are certainly compatible with the Bianchi identity V- Ein = and 
give rise, as we shall demonstrate below, to the expected Lorentz forces on 
charged matter due to vector fields. Since the Proca field couples to a current 
jq of Proca charged matter we have 

d^F+ (^^y ^a + ajq = 0. (5.5) 

The Proca charge current will be assumed to take the convective form 

jq = qng'kVg (5.6) 
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with constant Proca charge, q (of dimension length). 

In a similar manner we assume that the Proca neutral particle current is 
given by 

Jo = riQ-kVo. 

We postulate conservation of Proca charged particles 

d(nq -kV) = Q djg = (since q is constant). (5.7) 

If the neutral Proca matter is also conserved (as befits behaviour in the late 
post inflationary epoch) 

djo = (5.8) 

and one may then interpret the matter field equations as equations of motion 
for the fluid flows. 

To explicitly develop these equations we must compute the divergence of 
E. To facilitate this we introduce the tensors 

1 

= a ® a — - a{a) g (5.9) 



and 



so that 



= i,F F - ^ {F ^ ^ F)g (5.10) 



It is a standard result that 

V-^F^iuF (5.11) 

where the vector field U is defined by -kU = d-kF. To calculate the divergence 
of Eq consider 

V- (a a) = Vx„(q; ® Q;)(X") 

(5.12) 

o o o ^ ' 

= (V- a)a + a"- Vx^a = (*" d ^ a) a + a" Vx^a 
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where a" = a^X") and the relation V- /3 = ix-Vx,P = (-1)*'+^ *-^d^p for 
any p-form j3 has been used. Furthermore, for any 0-form / 



V-(/g) = d/ 

o 

since Vis metric compatible. Hence 

V- g) = d(a(a)) =e%Vx„(a(5)) 

= 2 (Vx.a)(a) = 2 a'^ Vx^a - 2 i~ d a 



(5.13) 



(5.14) 



where the vector field a is the metric dual of the form a. It follows from 
(13), (|5J2|) and (|J§ that 



V- Sc, = (^"M 7^ a) a + ij F. 
Together with ( p. Ill ) the divergence of S becomes 

C ?77,Q, \ ^ , _i , < 

- 1 (* d* a) a 



V-S = izF 



where the vector field Z is defined by 



*Z = d*F + 



•k a. 



Comparing ( |5.16|) and ( p3|) one has 

while the exterior derivative of both sides of ( [5. 51 ) yields 



h 



d-ka + crdj = 0. 



Therefore the divergence of S becomes 

o 

V- S = aqugivg F 



(5.15) 



(5.16) 



(5.17) 



using ( pJ[5D and (|5T^ ). 



For completeness we derive the divergence of the generic fluid stress: 
T ={c^p + P)V®V + Pg 
where g(V", V) = —1 and p = p{n, s). It follows from ( |5.13| ) that 

V-T= {c^ p + P)VvV + dP + V i.-^ d{{c^ p + P) i<V). (5.18) 

The third term on the right hand side of ( [5 .181 ) may be expanded as 

*-M((c2 p + P)^V)= ^-i{d(c2 p + p)^^V+(c^p + p)d^^V} 

= iv{d{c^P + P) + {c'p + P)i<~'di.V} ^' ' 

by using the relation 

for any 1-forms (3 and uj. Furthermore, the general conservation equation 
d(n 7t V") = dn A i<V + ndi<V = 

imphes that 

dikV = dnA kV = (iy dn) 1 

n n 

and so ^ 

i<~^di<V = (iy dn). 

n 

Thus ( ^.19| ) becomes 

^-M((c2 p + P)^V) = iv (d{c' P + P)- ^-^^^^ dnj . 

Applying the thermodynamic relation (|5.2|) to p{n, s) gives the divergence of 
T finally as: 

V- r = (c2p + ^) Vy\/ + nydP + iy(nTds)V^ 

where Ily dP = dP + (iydP)V^ is the transverse part of dP with respect to 
V. Thus the vanishing of the tangential components of the field equations 
and ( p. 41 ) gives rise to the isentropic conditions: 

Vo{so) = (5.20) 
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and 



V,{s,) = 



(5.21) 



while the vanishing of the transverse components yields the appropriate rel- 
ativistic Navier-Stoke's type equations for neutral and charged fluids respec- 
tively.: 

(c^ Po + Po) Vvo = -n^o dP (5.22) 



(c^ p, + P,) Vy, K = -liv, dP + 1^ iy, F. (5.23) 

In the case of negligible pressure {Pf/c^ ^ Pf^ / = 0; ^^e relation (|5.2| ) 
yields pj ~ A*/^/ where the constant p/ is identified as the Newtonian mass 
of the neutral and Proca charged particle respectively. In the limit where 
Einstein gravity is negligible we may introduce coordinates {t, x^, x^, x^} 
and write the metric g ~ — dt (S> dt + dx^ ® dx^ + dx^ ® dx^ + dx^ ® dx^. 
The Proca equation (|5.5| ) reduces to 



,4 



a%|; + |l_f£Ei)%.,,„-, (5.24) 



with the static ansatze ^) and a = cA{x^) dt for = 1, 2, 3. The 

Green function associated with ( ^.24|) is given by 



A = -- e-n^ (5.25) 

47rr 



where r = y (x^)^ + (x^)^ + (x^)^. It follows that a q A represents the Proca 
field produced by a single Proca charged source whose number density is 
approximated by the Dirac delta distribution nq(x^) = 6^{x''). If the world 
line of another particle has spatial velocity x'^(t) (x^/c ^ 1) then 

Vv',n«---|r+-"^^ + -»^- (5.26) 

If this particle also has Proca charge q and mass Pg then in the low energy 
and non-relativistic limit (|5.23|) and ( ^.25| ) imply that 



•■^ ac^q'^ d 
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Thus the "Yukawa potential energy" V{r) between these two Proca charged 
particles with separation r is 



Particles carrying the same sign of Proca charge attract or repel each other 
according as cr = +1 or — 1 respectively, with a force of range We give 
arguments below why it is of interest to consider the possibility that dark 
matter may interact with the Proca field having the polarity cr = +1. 



6. Implications from Galactic Dynamics 



The basic equations for our dark matter model are the Einstein equation 
(|5.1D, the matter field equations (|5.3|), (^^), the Proca equation (|5.5|), sup- 
plemented by the conservation equation (|5.7|) for the Proca current ( |5.6|) and 
the imposition of equations of state for thermodynamic variables satisfying 
the Gibb's relation ( ^.2| ). The dark matter is currently assumed to form 
an all pervading environment in which the less abundant ordinary matter 
moves. Since the ordinary matter is considered to have zero Proca charge its 
interaction with the dark matter is purely Newtonian in the non-relativistic 
low energy domain. We examine the extent to which this simple picture 
is sufficient to reconcile the galactic rotation curves with observation. For 
simplicity we shall treat the stars in a galaxy to be test particles that expe- 
rience Newtonian forces from a locally spherical distribution of low pressure 
[Pq <^ PqC^) dark matter gas at some temperature T and pressure Pq. 

Since the gas is charged one expects that its equation of state take the 
imperfect gas form ^2 



Pn 



nq + B2 nl + 0{nl) 



(6.1) 



where fc^ is Boltzmann's constant and B2 is the second virial coefficient given 
in terms of the Yukawa type potential ( ^.2?| ). 
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Under approximately adiabatic conditions ( |5.2| ) yields pq ~ PqUq with 
constant pq. The Newtonian dynamics of the self gravitating dark matter 
gas sphere gives rise to the Lane-Emden equation 

ATTr^(/^^= -AGpqM (6.2) 



where M = M(r) is the total mass within the sphere of radius r. If g 7^ 
then in terms of the dimensionless quantities b, r], u and N = N{u) defined 
by 

_ AV2nGpq 



(? q 



(6.3) 



and 



" = ^2 ^ (6-4) 

N = —M{r{u)) (6.5) 



6 = 32768^!^!i|!^ (6.6) 



?2f) becomes 



For M ^ 1 ( |6.71 ) has the asymptotic solution N{u) ~ corresponding 

to 



M{r)^'^^^r (6.8) 

irrespective of the sign of a. This gives rise to the asymptotic dark matter 
mass density distribution 

P«W = ^ 6.9 

2 7rG/igr^ 
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and the asymptotic dark matter pressure distribution 



27:0 fi^r^ 



(6.10) 



Thus 



^^3^ = jf^ and so the low pressure condition is satisfied if ksT -C 
fiqC^ (cold dark matter). Within such a dark matter gas, the rotation speed 
for stars, treated as test particles with zero Proca charge, following circular 
orbits of radius r under Newtonian gravity is approximately constant with 
value 



'GM(r) 



2kBT 



This relation may be used to eliminate 
dark matter density 



(6.11) 



op yielding the approximate 



47rGr2 



(6.12) 

220 Km/sec Equation ( |6l^ ) 



Taking a typical observational value 
gives a value Pq{r) ~ 6 x 10~^^ Kg/m^ of the dark matter density for r = 
lOkpc. 

Although there is no fundamental reason why the dark matter should 
have been in thermal equilibrium with the microwave background radiation 
at any time we may bound our estimates by assuming that its temperature 
now is no hotter than T = 2.7 K. For Vc ~ 220 Km/sec ( 6.11 ) implies that 

/ig ^ 1.5 X 10"22 Kg (6.13) 



or /i, ^ 10 G eV. 



7. Cosmological Models 

We now turn attention to the implications for cosmology by considering 
the class of homogeneous and isotropic Robertson- Walker type metrics 

g = dt®dt + (dx 0dx + f (dO dO + sin^ 6 d(p ® d0)) (7.1) 
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in local coordinates {t, Xi 0}; where S = S(t) denotes the scale factor and 
the function / = /(x) is specified by the curvature parameter k according to 



fix) = sinx? if k = 1, (spatially closed universe) 

fix) = Xj if k = 0, (spatially flat universe) 

fix) = sinhx, if k = —1, (spatially open universe). 

Both neutral and Proca charged fluids are supposed to have the common 
velocity: 

Vo = v, = i| 

c ot 

to a first approximation. Furthermore suppose that each fluid maintains the 
same entropy so that Sq and Sg are constants. This is consistent with the 
(weaker) isentropic conditions (|5.20| ) and (|5.21|) . The densities thus reduce 



to po = Poino) and pg = Pgirig) where uq = no(t) and Ug = Ugit). 
We adopt as ansatze for the Proca 1-form 

a = Adt (7.2) 
where A = Ait). It follows that F = da = and the Proca field equation 

d^F+(^^y ^a + aj^ = (7.3) 



gives 

a qugh^ 

i-t2i yy) 2 



A = ^-p^. (7.4) 



The Proca charge conservation equation ( p7^ ) yields the general solution: 

n, = ^ (7.5) 

where A^^ is a dimensionless constant. 

The Einstein equation ( |5.1| ) yields the differential equations 

)'-^(p-U^^ + c^A: = (7.6) 
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and ^ 

3|^ + 4(p + 3P-4Ca) = (7.7) 

where the total mass density and total pressure are defined by 

P = Po + Pq 

P = Po + Pq 



4_ 2^2 



and 

_ a Cm, 

Since ds/ = for / = 0, g, ( pl2|) becomes 

dp/ = /^/dn/ (7.8) 

where 

^ (7.9) 



It is instructive to write (p^) and ( |7. 7|) in terms of the Hubble parameter 

1 dSjt 
'S{t)~dt 



H(t).^?^ (7.10) 



and the deceleration parameter 



Then (p]^) and ( |7.7| ) become 



P = Cc. + ^-^ + ^^7^ (7.12) 
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Furthermore in terms of the critical density: 

Pc(t) = C + 1^ (7.14) 

we define 

^ - 4t = W- (7-15) 



Using ( [7. 121) one may ehminate (a in ( ^-ISD to get 



This parameter plays an important role in the discussion below. 



8. Dynamical Modifications to Friedmann Cosmologies 

To study the the effect of S on cosmology we first consider equations 
of state for both Proca charged and neutral fluids to be approximated by 
pressure-free dusts with conserved particle currents. It then follows from 
Po = 0, Pg = and (|7J) that 

Po = Po no 

Pq Pq 

with the constant chemical potentials fig, /io interpreted as the masses of the 
Proca charged and neutral particles. From the particle conservation equation 
hSD it follows that 

no = ^ (8.1) 



where A^o is a dimensionless constant. Hence from (|8.1| ) and (|7.5|) one has 



Po = ^ (8.2) 
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P, = ^ (8.3) 



oqhMq 



and ( [7.4|) becomes 

where Mq = Nq ^q, Mg = Ngfig. Using ( |7.8| ) it may be shown that ( |7.7| ) is 
satisfied provided S is not constant and the Equations ( [7.5|) and (|7.6| ) 

are satisfied. Now ( [7. 61 ) becomes 

^ Kr-7^ +c^ = (8-4) 



vat; 3 

where 

M = Mo + Mg. (8.5) 
In the following it is convenient to introduce the constant by 

e = -f-(ii^)\ (8.6) 

Note that the sign of is determined by the coupling polarity a. In terms 
of the dimensionless quantities defined by 

3c^t 

X 



SnGM 

27 c^e 



512 n^G^M^ 
and 

8nGM 

may be expressed as 

fdrV 1 1} 
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The solutions of this equation for various choices of (and k = 0) are illus- 
trated in Figures |l|, 0, |] around a; = where the solution with small is 
insensitive to the choice of k. 

For a = —1 where 6 < one finds solutions for Y exhibiting zeroes where 
the Levi-Civita geometry is singular as in the standard model of cosmology. 



For the case a = +1 where 9 > any positive valued solution to ( |8.71) has a 
non-zero minimum Y^i^ given by ^|y=y^.^ = 0. IfO<'(9<^l then 

~ (8.8) 

For Yyain non-zero the scale factor never vanishes and the geometry is non- 
singular. For example the Levi-Civita quadratic curvature invariant 

n = ^-'{R\^ ^R\) (8.9) 

has the value 

^_ 2A3c^ / 5 

~ ~5127c^G^M^ [TgY^ ~ Y^ ^ AY^ J ^' ^ 

which is clearly finite for y > 0. 

The invariant TZ has been scaled to 

- 512 .'G'M'^ 
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for display in Figures |I], 0, |^ in the vicinity of x = 0. We also show that in the 
domain Y(x) < Y{x) has an acceleration phase (^^ > , < 

obtained by differentiating ( |8.7| )). In Figure ^ the nature of the scale factor 
is exhibited over a larger region of x. The eternal nature of the cosmology 
is apparent for -(9 positive although depending on the size of S'min the model 
can readily accommodate hot dense phases where quantum conditions may 
be relevant. 
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9. Dynamical Modifications to Inflationary Friedmann 
Cosmologies 

In the above we have shown that the presence of the Proca field in the 
Einstein equations with a = 1 can modify the singular nature of the geometry 
associated with the Friedmann cosmologies. It is therefore of interest to 
examine how this field modifies these cosmologies when they are preceded by 
an infiationary expansion era. 

We recall the motivation for an infiationary phase starting at some t = ti 
after the standard singularity at t = 0. If tj is chosen to be the Planck time 
tpi then from ( [7.16|) a Friedmann radiation era predicts 



r ^ ^ ^ ^ ^ 10^° (9.1) 
Aj tpi 

where Aq = A(to) for to ~ 3.5 x 10^'^ sec and tpi ^ 5.3 x 10""'^ sec. It is 
accepted wisdom that unnatural constraints are required to maintain this 



ratio over such a long period. The standard infiationary model |34, Rq, 36 



|57| , PI attempts to bring the ratio r to order 1 by postulating a pre- 

Friedmann era from ti to tr in which the scale factor grows exponentially. 
Thus if the infiationary epoch ends at t = tj. with an e-folding e at an energy 
scale ~ 10^5 GeV: |ll 

_A,A,^^,„_, ^^^^ 

By choosing e > 60 the fine tuning puzzle is considered to be solved. 

We first note that the parameters of the purely infiationary models enable 
one to gain some initial estimates of G. Standard infiationary models pos- 
tulate that during the infiationary phase {t > ti) the Universe is dominated 
by the constant mass density p = p^, ~ 1 x 10^^ Kg/m^ which drives a rapid 



expansion of S until S* = 5,. ~ 10 m Thus for the pure infiation phase 
we prescribe the equation of state 

P=-c^p (9.3) 

so by differentiating ([T^) and using ( |7.7| ) the Proca modified Einstein equa- 
tions yields 
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and 

where the constant A is the "effective cosmological constant". (Standard 
inffation is described by (|9.5| ) with 6 = 0.) With ^ = at 5* = Smin then 

Q = ^(^Sl,-3kc'Sl\. (9.6) 



SttG 

Thus with log f #-) > 60 

\ '-'mill / 

S^i^ = S{U) < 10-2^ m (9.7) 

and ( p.6| ) imphes that 

e < 5 X 10"^^ Kgm^ (9.8) 

for /c = 0, 1, — 1. If the current Friedmann era of the Universe with t = 
to ~ 3.5 X 10^'' sec is dominated by dark matter with neghgible pressure 
{Pq{to) ^ Pqi'to) / c^) , then the current scale factor may be estimated to 
be S'(to) ~ 3 X 10^° light years = 3 x 10^^ m. Furthermore, 

Mg^M^ pM S{tof = 2 X 10^^ Kg (9.9) 

where the critical density has been chosen to be 6 x 10~^^Kg/m^. It then 
follows from (|J), (|03|) and 

Ng = ^ 1.3 X 10^^ m-2 (9.10) 



that 



2 / 1 \ 2 



16 G ( 1 e < lO-^*^^ mVKg2. (9.11) 



TUa'^ \hNq 



If ma is less than the Planck mass mpi = = 2.18x10 ^ Kg then |g| < 
2 X 10^^^ m which is much smaller than the Planck length ipi = J ^ = 1.6 x 
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10 111 and somewhat unnatural. Furthermore the dimensionless couphng 

"2 

^ associated with the effective coupling 

is many orders of magnitude smaller than the fine structure constant. The 
origin of these values is the high e-folding e needed in the standard inflation- 
ary model to accommodate the flatness problem associated with the standard 
Friedmann curvature singularity. However in the model with Proca fields 
there need be no such singularity and one may adjust the size of the mini- 
mum scale factor to accommodate the flatness problem with a smaller e. We 
estimate the necessary parameters by requiring that q ~ £pi and nia ~ m^\. 
In this case the dimensionless coupling strength between two Proca charges 
is 

e 1 1 



he 327r2 315 



(9.13) 



which is of the order of fine structure constant. 

For simplicity we arrange that at the end of the inflationary epoch the 
Universe proceeds to the current time in a Friedmann era dominated by dark 
matter modelled as a Proca charged fluid with negligible pressure. From ( |9.6| ) 
the minimal value of the scale factor S'min = 'S'(O) may be determined by 
and A 



5mi„ = V2 [-^) (9.14) 

neglecting k. If the inflationary era ends at time U with a scale factor Sr = 
S{tr) and 

Pinflation(^r) PFriedmann {tr) (9.15) 

one has 

A = (9.16) 

Or 



With the value of Nq from ( |9.1CI|) and g ~ i'pi , ma ~ mpi, (|8.6|) gives 

e = 3xl02^Kgml (9.17) 
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It follows from {^A§, and ( pJ[5| ) that 



r=^^ = ^(^y (9.18) 



Eliminating ^min using (|9.14|) and assuming the desired value r < 1 we may 



solve for 5,: 

Sr>2x 10^ m (9.19) 
Inserting this value in ( |9.14|) yields the value of the minimal scale factor 

5mm > 10 cm (9.20) 

This corresponds to e = log {Sr/ •S'min) > 20 which is a considerably smaller 
bound than the traditional inflationary expansion based on Einsteinian grav- 
itation alone. 



10. Conclusion 

Motivated by the structure of a class of actions that involve (in addi- 
tion to the generalised Einstein-Hilbert action) terms including the torsion 
and metric gradient of a general connection on the bundle of linear frames 
over spacetime, the consequences of Einstein-Pro ca gravitation coupled to 
matter have been examined. This theory is written entirely in terms of the 
traditional torsion free, metric compatible connection where all the effects of 
torsion and non-metricity reside in a single vector field satisfying the Proca 
equation. In such a theory the weak field limit admits both massless tensor 
gravitational quanta (traditional gravitons) and massive vector gravitational 
quanta. The mass of the Proca field is determined by the coupling con- 
stants in the parent non-Riemannian action. The interaction mediated by 
the new Proca component of gravitation is expected to modify the traditional 
gravitational interaction on small scales. In order to confront this expected 
modification with observation we have constructed an Einstein-Pro ca-Fluid 
model in which the matter is regarded as a perfect thermodynamic fluid. 
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We have suggested that in addition to ordinary matter that couples gravita- 
tionally through its mass the conjectured dark matter in the Universe may 
couple gravitationally through both its mass and a new kind of gravitational 
charge. The latter coupling is analogous to the coupling of electric charge 
to the photon where the analogue of the Maxwell field is the Proca field 
strength (the curl of the Proca field). If one assumes that the amount of 
dark matter dominates over the ordinary matter in the later phase of evolu- 
tion of the Universe, that the Proca field mass is of the order of the Planck 
mass and the appropriate coupling to the dark matter is of the same or- 
der as the fine structure constant then one finds that such hypotheses are 
consistent with both the inflationary scenario of modern cosmology as well 
as the observed galactic rotation curves according to Newtonian dynamics. 
The latter follows by assuming that stars, composed of ordinary (as opposed 
to dark matter), interact via Newtonian forces to an all pervading back- 
ground of massive gravitationally charged cold dark matter in addition to 
ordinary matter. The novel gravitational interactions are predicted to have 
a significant infiuence on pre-infiationary cosmology. For attractive forces be- 
tween dark matter charges of like polarity the Einstein- Proca-matter system 
exhibits homogeneous isotropic eternal cosmologies that are free of cosmolog- 
ical curvature singularities thus eliminating the horizon problem associated 
with the standard big-bang scenario. Such solutions do however admit dense 
hot pre-infiationary epochs each with a characteristic scale factor that may 
be correlated with the dark matter density in the current era of expansion. 

The Einstein-Proca-Fluid model offers a simple phenomenological de- 
scription of dark matter gravitational interactions. It has its origins in a 
geometrical description of gravitation and the theory benefits from a varia- 
tional formulation in which the connection is a bona fide dynamical variable 
along with the metric. The simplicity of the model is a consequence of the 
structure of a class of non-Riemannian actions whose dynamical consequences 
imply that the new physics resides in a component of gravitation mediated 
by a Proca field. It will be of interest to confront the theory with other 
aspects of astrophysics such as localised gravitational collapse, the nature of 
the inflation mechanism and the origin of dark matter. 
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Figure 1: The Early Universe for > 



The solution to (|8.7|) with a small positive value of ■& has an asymptotic 



form Y{x) ~ d^^^ + l^?"^/^ around x = 0. Such a Y{x) has a minimum 
at a; = 0, where the corresponding dimensionless curvature invariant ^(s) 
defined in ( ^.11| ) has a maximal value. Note that Y{x) has an acceleration 
phase {Y{x) > 0) within the region -2 (^9/3)^/2 < a; < 2 (^9/3)^/^ 
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Figure 2: The Early Universe for -(9 = 
The solution Y{x) ~ |3x/2p/^ for the standard Friedmann model is obtained 
by choosing 'd = in (|8.7] ). The "big-bang" or "big-crunch" singularity 
follows from Y ^ and S — > oo as x — >■ 0. Y{x) is negative definite 
throughout spacetime. 
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Figure 3: The Early Universe for < 
With a choice of negative ( |8.7] ) yields the asymptotic solution Y{x) ~ 
13x1^/^ which approaches zero as x — more rapidly than Y{x) de- 
scribed by the Friedmann solutions, thus leading to a more severe singular 
behaviour of the spacetime curvature. 
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Figure 4: Eternal and Oscillating Cosmologies 
The global behaviour of Y{x) for i) = 0.1 exhibits three types of eternal cos- 
mology. If the curvature parameter A; = 1 an oscillating universe is depicted 
by a periodic Y(x). If /c = or —1 the universe contracts to a minimal scale 
and then expands to infinity. In all three cases the Levi-Civita curvature is 
regular throughout spacetime. 
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